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Abstract. In quantum cosmological models, constructed in the framework of Friedmann-Robertson- 
Walker metrics, a nucleation of the Universe with its further expansion is described as a tunneling transition 
through an effective barrier between regions with small and large values of the scale factor a at non-zero (or 
zero) energy. The approach for describing this tunneling consists of constructing a wave function satisfying 
an appropriate boundary condition. There are various ways for defining the boundary condition that lead 
to different estimates of the barrier penetrability and the tunneling time. 

In order to describe the escape from the tunneling region as accurately as possible and to construct the 
total wave function on the basis of its two partial solutions unambiguously, we use the tunneling boundary 
condition that the total wave function must represent only the outgoing wave at the point of escape from 
the barrier, where the following definition for the wave is introduced: the wave is represented by the wave 
function whose modulus changes minimally under a variation of the scale factor a. We construct a new 
method for a direct non-semiclassical calculation of the total stationary wave function of the Universe, ana- 
lyze the behavior of this wave function in the tunneling region, near the escape point and in the asymptotic 
region, and estimate the barrier penetrability. We observe oscillations of modulus of wave function in the 
external region starting from the turning point which decrease with increasing of a and which are not shown 
in semiclassical calculations. The period of such an oscillation decreases uniformly with increasing a and 
can be used as a fully quantum dynamical characteristic of the expansion of the Universe. 
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1 Introduction 

Tracking the evolution of the Universe back in time on the basis of the classical equations of general relativity, one 
encounters singularities, where such equations break down. In order to understand what really happens in the formation 
of the Universe, many people came to the point of view that a quantum consideration of this process is the deeper one. 
So, according to [1], in the quantum approach we have the following picture of the Universe creation: a closed Universe 
with a small size is formed from "nothing" (vacuum), where by the word "nothing" one refers to a quantum state 
without classical space and time. Here, the wave function can be used for a probabilistic description of the creation of 
the Universe and its subsequent expansion. The first papers with the quantum approach for the description of Universe 
formation and its initial expansion may be [SIE], and shortly afterwards many other papers appeared in this field, 
pointing to a rapid development of the quantum approach in cosmology (for example, see H , 5 , 6 , 7 , 8Tf^ [TUlfTTlfT^ and 
some discussions in [13 1 114 ] and references therein). 

A key point in the construction of the proper wave function is the choice of a boundary condition. We find a 
natural and clear condition in the papers of Vilenkin — the tunneling boundary condition pilTB] , according to which 
the wave function in the region of the escape point from the barrier must represent only an outgoing wave. The 
simplest answer to this question is given in the semiclassical approach (for example, see Refs. [13pi4j). However, the 
validity of the semiclassical approximation decreases when we pass closer to the escape point, which is exactly the 
point where we want to impose the boundary condition, and we must take into account exact solutions here. It is 
natural to suppose that the direct quantum approach is richer, accurate and detailed. And perhaps, it should be more 
effective going beyond the semiclassical study of nucleation of the Universe and its further evolution. On the other side, 
the semiclassical methods are practically shown as accurate enough in a quantum calculation of the rates and other 
characteristics of the evolution of the Universe in cosmological models. The semiclassical methods are more stable and 
convergent in computer calculations (which is very important for obtaining reliable values of the calculated parameters 
in cosmological scales). At present, a semiclassical power formalism has been developed where the physical meaning of 
the parameters and characteristics is clear enough intuitively [T5] . On such a basis, one can ask the question: whether 
is the direct non- semiclassical quantum approach for the determination of the wave function of the Universe really 
better in comparison with semiclassical methods? 

Let us find such characteristics, which are less effectively determined by the semiclassical approach. In order to do 
this, let us consider the rules of correspondence between two semiclassical wave functions in the tunneling region and 
in the above barrier region close to the turning point up to the second correction of the semiclassical approximation 
(for example, we use (47,5) in p. 208 and (50,2) in p. 221 of "16]). We conclude to the following. 

— The semiclassical wave function thus defined in the tunneling region is real and we obtain a zero flux in this region. 
In other words, we have no any propagation of waves under the barrier. So, the semiclassical approximation is 
less effective for study of tunneling processes with non-zero fluxes. But this is a case of decay of the quantum 
system (contrary to elastic scattering in nuclear theory) and a similar process is used in the quantum description 
of the nucleation of the Universe and its evolution at its first stage. The third correction of the semiclassical 
approximation (applied in a more restricted spatial region) transforms a real wave function into a complex one 
(see, for example, (46,11) in [16j) and now we see that such a picture without propagation of the wave under 
the barrier is only approximation. So, the first conclusion is: the semiclassical approaches are less effective for a 
(detailed) study of tunneling processes with strong propagation of waves (i. e. for such spatial regions under barrier 
where the phenomenon of tunneling is present "strongly" ) . 

— The semiclassical wave function in the above barrier region is complex. From this we obtain non- conservation 
of the flux in the transition through the boundary at the turning point. Here, the flux calculated in the second 
semiclassical approximation describes propagation of the Universe in a "classical sense" but starting from the 
turning point, while exact calculations give conservation of the flux in the whole spatial region starting from zero. 

— The flux non-conservation breaks down the correct change of the phase of wave function along the spatial axes 
(one may suppose that this violation can be essential). So, we come to the next characteristic — the phase of the 
wave function — which can be less correctly estimated in the semiclassical approximation. 

— From this we come to the phase time, which is deflned on the basis of the phase of the wave function and which 
is used in quantum dynamical theories. Thus, we come to the supposition that the semiclassical approach is less 
effective in a study of quantum dynamics (especially in the tunneling region). 

— The phase includes information about interaction between barrier and wave, and it can be interesting in a detailed 
study. 

— The boundary condition is imposed in the non-semiclassical region while we construct the outgoing wave in the 
semiclassical limit. One can suppose that the connection between the boundary condition and the formulation of 
the wave (for the same exact potential form) should be realized more accurately in the direct quantum approach. 

Therefore, in order to estimate the probability of the nucleation of the Universe in the framework of the quantum 
approach with the highest accuracy, the necessity of the construction of a direct (non-semiclassical) method for an 
unambiguous construction of the wave in the region close to the escape point has become clear. A possible answer to 
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the question above may be by the semiclassical approach some characteristics are determined very well, other ones 
can be determined less accurately and correctly. 

In this paper new strict and weak definitions of the wave propagating inside a potential of arbitrary shape and 
construction on this basis of the total wave function for cosmological potentials are proposed. After a short description 
of a simple variant of the quantum model of the evolution of the Universe in Sec. a new method of the direct 
(non-semiclassical) determination of the stationary wave function at zero energy is presented and its behavior in the 
barrier region, near the escape point and in the asymptotic region, is analyzed in Sec. [3l and the barrier penetrability 
by the proposed approach is estimated with comparison with semiclassical calculations in Sec. |31 



2 Cosmological model in the Friedmann-Robertson-Walker metric 

One can describe the nucleation of the Universe by a simple model. Here, we shall consider the simplest variant of the 
homogeneous and isotropic Universe. We consider the Friedmann-Robertson-Walker (FRW) metric which defines the 
most general form of 4-dimensional spacetime with spherical spatial symmetry (see Ref. p. 438; also see Refs. [HI 

ds^ = ~dt^ + a^{t) ■(^+ r^{de^ + sin^ 9 dc^^)] , h{r) - 1 - kr^ , (1) 
\h{r) J 

where t and r, 9, (j> are time and space spherical coordinates, the signature of the metric is +, +, +) as in Ref. [20] 
(see p. 4), a{t) is an unknown function of time and fc is a constant, the value of which equals +1, or —1, with 
appropriate choice of units for r. Further, we shall use the following system of units: h ^ c = 1. For k = —1, the 
space is infinite (Universe of open type), and for k = +1 the space is finite (the Universe of closed type). For fc = 1 
one can describe the space as a sphere with radius a(t) embedded in a 4-dimensional Euclidian space, and the function 
a{t) is referred to as the "radius of the Universe" . In the general case, it is called the cosmic scale factor. This function 
contains information of the dynamics of the expansion of the Universe, and therefore its determination is an actual 
task. 

One can find the function a{t) using the Einstein equations in this metric (we use the signs according to the chosen 
signature, as in Ref. |20] p. 8; the Greek symbols fj. and denote any of the four coordinates t, r, 9 and 0): 

R^..-^g^..B = 8^TGT^,, (2) 

where i?^^ is the Ricci tensor, R is the scalar curvature, T^^ is the energy-momentum tensor, and G is Newton's 
constant. 

Substituting the components of the Ricci tensor the scalar curvature R, and the components of the energy- 
momentum tensor T^i, (see Ref. |20j . p. 8) into the equation ([2]) at /i = i/ — 0, we obtain the Friedmann equation (see 
p. 8 in Ref. [20]; p. 3 in Ref. [H]; p. 2 in Ref. T]): 



=■2 ' fc - H'^a^ = 0, (3) 



a 



where 



H=\^-P, (4) 



and p is the energy density. Equation ([3]) determines the function a(t), describing the classical dynamics of the extension 
of the Universe, as determined by the energy density p. 

Let us define the action for the model. Usually, in the construction of the cosmological model the action must 
include both the geometry (curvature) of spacetime and matter fields. There are various papers proposing different 
variants of the cosmological models with inclusion of different types of the matter fields. In the given paper, in the 
construction of the action we restrict ourselves to the simplest component of the matter fields — the vacuum only, 
directing the main attention to the development of the method for calculation of the wave function. We define the 
action as in Ref. [T] (see (1), p. 2): 

R \ , 4 „ 6a^ + Gad + 6fc 



Substituting the scalar curvature, rewriting the energy density p through H according to and integrating the d 
term in the resulting expression by parts with respect to t, we obtain the lagrangian (see Ref. [I], (11), p. 4): 

L{a,d) = ^(^k~d^ -H^a^y (6) 
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Here, we shall consider the variables a and d as generalized coordinate and velocity, respectively. Let us find the 
generalized momentum conjugate to a: 

dL{a, a) . 

Pa = TT. = -aa- (I) 

da 

Further, we shall study the Universe of closed type only (k — +1). We then obtain the hamiltonian: 

h{a,pa) ^pa~ L{a,a) ^ -^(pl + - H^a'^y (8) 

This hamiltonian describes the classical dynamics of the Universe, which serves as a basis for classical cosmological 
models. 

The passage to the quantum description of the evolution of the Universe is obtained by the standard procedure of 
canonical quantization in the Dirac formalism for systems with constraints. We do the following. 

— We suppose that the variables a, pa and functions of them (for example, the hamiltonian) are operators, acting on 
a new function of state, named the wave function of the Universe, ip{a). 

— We change the generalized coordinate a and momentum pa into operators: 

a ^ a ^ a, Pa ^ Pa ^ -i-^. (9) 

oa 

Here, the rule of the canonical correspondence between classical and quantum Poisson brackets is fulfilled. 

— In order to take into account fS]) in the construction of the quantum model, we consider it as a constraint imposed 
on this quantum system, which gives us the eigenvalues (and eigenfunctions) of the operators which correspond 
most closely to fulfillment of this equation in the quantum case. 

As a result of the quantization, we obtain the Wheeler-De Witt (WDW) equation (see Ref. [1], (16)-(17), in p. 4, 
Ilir22j): 

'^+a'-HV^ (p(a)=0. (10) 

One can see that this equation looks similar to the one-dimensional stationary Schrodinger equation on a semiaxis (of 
the variable a) at zero energy with potential 

V{a)=a^-H^a'^. (11) 

Taking into account that in different papers different coefficients of and are used in the potential (fTTj) (as a result 
of a different construction of the action), for convenience, we rewrite this potential in a generalized form: 

V{a) = Aa^ ~ Ba^. (12) 

So, in particular, for the choice of parameters in Ref. |23| 

A = 36, S = 12 • yl, yl = 0.01, (13) 

the potential is shown in Fig. [TJ Let us denote the point of intersection between the potential V{a) and the a axis as 




Fig. 1. The potential V{a) with parameters A = 36, B = 12 ■ yl and A = 0.01 (point atp = 17.320508). 
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atp. Since we shall be studying the process of the escape from the barrier at zero energy, let us call the point atp the 
escape point. 

After the passage from the classical description of the Universe evolution to the quantum one, all the information 
of the dynamics of the Universe expansion in the early stage and its further evolution, which was previously contained 
in the lagrangian (or hamiltonian), passes into the wave function <p(a). Therefore, many papers today are devoted to 
the study of this function and any approach to its determination with the highest accuracy may be of interest. 



3 Determination of the wave function 

Today, the two most prevailing ways for the determination of the wave function ip{a) are the Feynman formalism of path 
integrals in multidimensional spacetime, developed by the Cambridge group and other researchers, called the "Hartle- 
Hawking method" ^), and a method based on direct consideration of tunneling in 4-dimensional Euclidian spacetime, 
called the "Vilenkin method" [4ll9l [T0lfT3] ) . We shall use the second way, which looks more natural and clearer for us. 
Here, the prevailing approach for the description of tunneling consists in its semiclassical consideration (for example, 
see Ref. [13ll4,24,25 ). An attractive side of this approach is its simplicity in the construction of decreasing and 
increasing partial solutions for the wave function in the tunneling region, the outgoing wave function in the external 
region, and the possibility to define and estimate in a simply enough way the penetrability of the barrier, which can 
be used for obtaining the duration of the nucleation of the Universe. We shall be interested in another way of direct 
non-semiclassical determination of the wave function (/3(a), in order to describe most accurately the process of the 
escape from the barrier, which describes the formation of the Universe and its further expansion at the first stage. 

3.1 The form of the wave function close to the escape point 

For finding the wave function ip{a) close to the escape point atp, we expand the potential V{a) in ()12p in powers 
q = a — atp (for small q) and restrict ourselves to the linear term only: 

Viq) = Vo + Viq. (14) 

From (fT2l) we have 



Vq — y(a = atp) = Aa 



Vi = 2atp {A - 2Ba%) = 2J - ■ {A - 2A) 





(16) 



and the Schrodinger equation with the potential V{q) with the new variable q becomes 

_d^ 
dq 



— ^{q)^\V,\q^{q) = 0. (17) 



Changing the variable to 
(fT7|) is transformed into 



i = \Vi\'/\, (18) 

^^iO+^fiO^O- (19) 

From quantum mechanics we know two linearly independent exact solutions for the function ip{^) in this equation 
— these are the Airy functions Ai (^) and Bi (^) (for example, see Ref. [321 , p. 264-272, 291-294). Expansions of these 
functions into power series at small ^, their asymptotic expansions at large |^|, their representations through Bessel 
functions, zeroes and their asymptotic expansions are known. We have some integrals of these functions, and also the 
form of the Airy functions in the semiclassical approximation (which can be applied at large |^|). In some problems of 
the analysis of finite solutions ip{S,) in the whole range of ^ it is convenient to use the integral representations of the 
Airy functions (see (10.4.32) in Ref. [26], p. 265, a = 1/3, taking into account (10.4.1)): 

+ 00 



Ai (±0 = ^ y cos T du, 



Bi (±0 = ^ / y S.y^ + sin T ^uj 



(20) 



du. 
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Furthermore, we shall be interested in the solution of the function (p{^) which most accurately describes the outgoing 
wave in the range of a close to the point atp. However, it is not clear what the wave looks like near the point atp in 
the potential studied in general, and which linear combination of the functions Ai (^) and Bi (^) defines it in the most 
accurate way. 

The clearest and most natural understanding of the outgoing wave is given by the semiclassical consideration of 
the tunneling process. However, at the given potential the semiclassical approach allows us to define the outgoing 
wave in the asymptotic region only (while we can join solutions in the proximity of atp by Airy functions). But it 
is not clear whether the wave, defined in the asymptotic region, remains outgoing near the atp. During the whole 
path of its propagation outside the barrier the wave interacts with the potential, which must inevitably lead to 
a deformation of its shape (like to appearance of a phase shift in the scattering of a wave by a radial potcmtial 
caused by interaction in scattering theory). Could it turn out that the potentials used in cosmologic:al models give 
a significantly larger change of the shape of the wave caused by interaction in a comparison with the potentials used, 
for example, for the description of nuclear collisions in the framework of scattering theory? Moreover, for the given 
potential there is a problem with obtaining convergence in the calculation of the partial solutions for the wave function 
in the asymptotic region. According to our calculations, a small change in the range of the definition of the wave in 
the asymptotic region leads to a significant increase of errors, which requires one to increase the accuracy of the 
calculations. Therefore, we shall be looking for a way of defining the outgoing wave not in the asymptotic region, 
but in the closest vicinity of the point of escape, atp- In a search of solutions close to the point atp, i. e. at small 
enough |^|, the validity of the semiclassical method breaks down as |^| approaches zero. Therefore, we shall not use 
the semiclassical approach in this paper. 

Assuming the potential V{a) to have an arbitrary form, we define the wave at the point atp in the following way. 

Definition 1 (strict definition of wave) The wa,ve is such a, linear combination of two partial solutions of the wave 
function that the change of the modulus p of this wave function is closest to constant under variation of a: 



'2 



0. (21) 



According to this definition, the real and imaginary parts of the total wave function have the mutually closest behaviors 
under the same variation of a, and the difference between possible maximums and minimums of the modulus of the 
total wave function is the smallest. For some types of potentials (in particular, for a rectangular barrier) it is more 

convenient to define the wave less strongly. 

Definition 2 (weak definition of wave) The wave is such a linear combination of two partial solutions of wave 
function that the modulus p changes minimally under variation of a: 



0. (22) 



According to this definition, the change of the wave function caused by variation of a is characterized mainly by its 
phase (which can characterize the interaction between the wave and the potential). 

Subject to this requirement, we shall look for the solution for the function (p{^) in the following form: 

Umax 

m= J expz(-y+/(Oii) du, (23) 



where A/' is a normalization factor, f{£_) is an unknown continuous function satisfying f{£_) const as ^ ^ 0, and Umax 
is the unknown upper limit of integration. In such a solution, the real part of the function /(^) gives a contribution 
to the phase of the integrand function only, while the imaginary part of /(^) deforms its modulus. 

Let us find the first and second derivatives of the function ^{^) (a prime denotes a derivative with respect to ^): 



"max g 

^^iO=i I f'uexpi(-^ + f{Ou)du, 



- / {^f"u - iffu') cxp^(-y + fiOu) du. 



(24) 
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From this we obtain: 



du. 



(25) 



Considering solutions at small enough values of |^|, we represent /(^) in the form of a power series: 



(26) 



Tl = 



where /„ are constant coefficients. The first and second derivatives of /(^) are 



4-00 



(27) 



Substituting these solutions into (|25l) . we obtain 

d2 



2'Z'(e)+C'z^(C) 







(2iu /2 - ^2 /2) + (6m h - ^u"" fj2 + l) e+ 



+ CXJ r n 

+ J2 iu {n + l){n + 2) /„+2 - u'^ J2 i'^ - m + l)(m + 1) /„-„i+i/m+i 

n=2'- m=0 

,3 



(28) 



X exp I ^ 

Considering this expression at small |^|, at the first step we use the following approximation: 



expj(-y + fuj expi(-y + fgl 



(29) 



Requiring the condition (|28p to be satisfied for different ^ with different powers n, we obtain the following system: 
e ■■ J (2m /2 - f!) exp * (- y + /o") du = 0, 



: y" (6m h - /1/2 + l) exp^(-y + /o") = 0, 



(30) 



m (n + l)(n + 2) /„+2 - E ~ + + 1) fn-m+lfm+l 



m=0 



X exp i — + /qw j — 0. 
Assuming the coefficients /o and /i to be given, we find the following solutions for the unknown /2, /s and /„: 



f? J2 



h 



4/1/2 J2 - Jo 



2j Ji ' 6i Ji 

n 

{n~m + \){m + 1) fn-m+lfm+l 



m— 



/n+2 



i (n + l)(n + 2) 



J2 



(31) 



(32) 
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where the following notations for the integrals have been introduced: 



Jo = / cxpi{ — j + fou] du, 



Ji= I uexpi[ — — + faujdu, (33) 








expi(^ — — + fouj du. 



Thus, we see that the solution for the function cp{£^) has arbitrariness in the choice of the unknown coefficients 
/o, /i and the upper limit of integration Umax- However, the solutions found, PT|) and ([5^ . define the function /(^) 
so as to ensure that the equality (|19p is exactly satisfied in the region of a close to the escape point atp. This proves 
that the function (/?(^) in the form \2SI\) with an arbitrary choice of fo, fi and Mmax is the solution of the Schrddinger 
equation i23\) near the escape point atp and its accuracy is maximal at a = atp. In order to bring the solution •f'(C) 
into the well-known form of the Airy functions, Ai (^) and Bi (^), we select 

fo - 0, /i = 1. (34) 

At such choice of the coefficients fa and /i , the integrand function in the solution ([23|l up to has a constant modulus 
and a varying phase (the coefficient /2 deforms the modulus, but it is fulfilled at ^•^). Therefore, one can expect that 
the solution ([23]) at the escape point atp describes the wave with the proper shape. 



3.2 Calculation of two partial solutions of the wave function 

Not knowing the behavior of the wave function and its derivative as a function of a, we shall be looking for their two 
partial solutions independently in the following way. At first, we shall define the wave function and its derivative at 
a selected point, and we shall calculate the wave function and its derivative in the region close enough to this point 
using the method of starting of the solution; this is presented in the next two sections. Here, for the partial solution, 
which increases in the barrier region, we use the starting point to be a = 0, and for the second partial solution, which 
decreases in the barrier region, we select the starting point to be the escape point atp. Further, we shall calculate the 
wave function and its derivative in the whole required range of a using the method of continuation of the solution, for 
which we select the Numerov method with a constant step. 



3.2.1 Determination of the wave function close to zero. A regular solution 



We shall be looking for the regular solution for the wave function close to a = 0. Let's write the wave function close 
to a = in the following form: 

(/3(a) = ci ^ a„ a", (35) 

where a„ are constant coefficients. 

One can find the unknown a„ from the Schrodinger equation with the potential (|12p . From psp we obtain the first 
and second derivatives of the wave function: 

ip'{a) = ci X) On = ci + 1) «n+i a") 

n=l n=0 

ip"{a) =ciJ2 (" + l)na„+i a""^ ^ ci J2 (n + 2)(n + 1) a„+2 a". 

n—l n— 

The Schrodinger equation and the potential have the forms 

-if" (a) + V{a) ip{a) = 0, V{a) =Aa^-Ba'^. (37) 
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Substituting (|35p and (|36p for the wave function and its second derivative, we obtain 



or 



+00 

V"(a) = ci ^ (n + 2) {n + 1) a„+2 a" = 

ri=0 



+ 00 



(p{a) = [Aa^ — B a^j ci J2 anO^ ^ 

-f 00 -\-OQ 



+ cxj +00 

ci -j A ari-2 a" - B a„_4 a" 

n=2 n=4 



Ci< ^^ao a 



2 -L fli 



+00 , 

+ [n + 1) an+2 = A ia^ a? + ai of 
n=o ^ 



+ CXJ 



Let us write the expressions at a" with the same powers n: 



n 


= 


-> 2 • 


1 


• 02 • a° 


= 0; 




n 


= 1 


^ 3 • 


2 


• 03 • fli 


= 0; 




n 


= 2 


-> 4- 


3 


• 04 • a'^ 


= A • ao 


a'; 


n 


= 3 


^ 5 • 


4 


• 05 • a'^ 


= A • ai 




n 


> 4 


^ (n 


+ 


2) (n + 


1) • an+2 


• a" = (^a 



(38) 



(39) 



Thus we obtain recurrent relations for the calculation of the unknown a„ 



02 



0, as = 0, a4 



Aap 
12 ^ 



05 



Aai 



■, an+2 



^an-2 - BUn- 



(n + l)(n + 2) 



at n > 4. 



Given values for ao and ai, using (I40p one can calculate all a„ needed. 
Analyzing ((35)) at zero, we find (here, we use ci = 1): 



(40) 



(41) 



So the coefficients oo and oi determine the wave function and its derivative at zero. Let us write all possible cases: 

(^(0) = 0, (^'(0)^0^00 = 0,01^0, 

(^(0)^0, (^'(0) = -> 00 7^0,01=0, (42) 
^{0) + 0, (^'(0) ^ ^ 00 ^ 0, 01 ^ 0. 

One can consider expressions (j42|) as different boundary conditions for the wave function and its derivative at zero. 
Implying two different boundary conditions through oo and oi (in such a way that they locate the first node for the 
wave function at different places), we shall obtain two linearly independent partial solutions >pi{a) and </32(a) for the 
wave function close to zero. At the present stage, not knowing the coordinates of the maxima and coordinates of the 
next nodes for the partial solutions for the wave function, we select 



^i(O) = 0, ^'i(O) ^ ^ oo = 0, 01 = 1, 
¥^2(0) 7^0, ^^(0) = 0^00 = 1,01=0. 



(43) 



However, analysis has shown that the error in the calculation of the decreasing component of the wave function can 
grow significantly when a increases (depending on the choice of the parameters of the potential), which makes the 
calculation of that component inefficient. Therefore, in this way we shall obtain one partial solution for the wave 
function only, which increases in the barrier region (we select the first condition from (j43p ). 



3.2.2 Determination of the wave function near an arbitrary selected point 

Now we shall find a regular solution for the wave function near a point a^, which can be selected arbitrarily on the 
semiaxis. We write such a wave function in the following form: 

+00 +00 
(/j(a) = C2^5„ (a- o^)" = C2^6„a", (44) 

n=0 n=0 
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where 6„ are new constant coefficients, and we introduced the new variable 



a — a ~ ax- 

Let us find the first and second derivatives of the wave function ((44|) : 



+00 



(p'{a) = C2 X] " a" ^ = C2 X] (" + 1) 



71=1 

+00 



(p"ia) = C2 X ("+ ^)nb„+i a 



We rewrite the potential through the variable a: 



n=0 
n-1 



+00 

= C2 E(« + 2)(n + l)6„+2 a". 

n=0 



(45) 



(46) 



where 



V{a) ^ Aa^ - B a-^ ^ a'^{A - B a^) ^ {a + a^f^A -B{a + a^f) 
= {al + 2a^a + a2)(A - B al - 2B a^a - Bo?) = 
= Co + Ci a + C2 52 + C3 + Ci a^, 

Co = Aal-Bal, 

Ci = 2ax{A- Bal) -2Bal^ 2Aax-4:Bal, 
C2 = A - B al - iB al - B al = A - 6B al, 
C3 = -2Bax -2Bax = -4:Bax, 
C4 = -B. 



(47) 



(48) 



Substituting the wave function (|44|), the second derivative of the wave function ([46]) and the potential ([47]) into 
the Schrodinger equation, we obtain 

+00 



V"{a)^c2 E(« + 2) (n+l)6„+2 a" 



n=0 
.2 u „4 



+00 



= (Aa^-Ba^) ip{a) = Cq + Ci a + C2 + C3 + C4 a'' C2 E = 
r +00 +00 +00 +00 +00 

= C2 Co E &n a" + Ci E ^« 5"+^ + C2 E fen + C3 E fen + C4 E fen 

L n^O n^O n^O n^O n^O 

r +00 +c)o +00 4-CX3 +00 

= C2 Co E fen a" + Cl E fen-l + C2 E fen-2 + C3 E fen-3 + C4 E fen-4 5" 
L n—O n—1 n—2 n— 3 n— 4 



C2 |Co 5o + (Co 61 + Cl 5o) a + (Cq 62 + Ci 61 + C2 60) + (Co 63 + Ci 62 + C2 bi + C3 5o) a^+ 

+ 00 

- E (Co fen + Cl 5„_i + C2 bn-2 + C3 5„-3 + C4 6„-4) 
n=4 



or 



+00 

E (« + 2) (71 + 1) 6„+2 a" = 

= Co 5o + (Co 61 + Cl 5o) a + (Cq 62 + Ci 5i + C2 60) a^ + (Co 63 + Ci 62 + C2 bi + C3 bo) 0?+ 

+ 00 

+ E (Co fen + Cl bn-l + C2 5„-2 + C3 6„-3 + C4 6„-4) S". 

n=4 



(49) 



We write the expressions at a" with the same powers n: 



n 


= 


^ 2 • 


1 


•62-50 


= Co6o 


n 


= 1 


^ 3 • 


2 


•63-51 




;Co6 


n 


= 2 


^ 4- 


3 


• 64 • 




Co 6 


n 


= 3 


-> 5 • 


4 


• 65 • 53 




;Co6 


n 


> 4 


->(n 


+ 


2) (n + 


1) 


■ 6„+ 



(Co 63 + Cl 62 + C2 61 + C3 60) • a^; 

■ bn+2 ■ a" = (Co 6„ + Cl 6„_i + C2 6„-2 + C3 6„_3 + C4 6„_4) • a". 
As a result, we obtain recurrent relations for the determination of the unknown 6„: 

, C060 , C061+C160 , C062 + C161 +C260 , C063 + C162 + C261 + C360 

62 = —7;—, 63 = 7i ' fe4 = 77^ 7 fes = 



(50) 



'J71+2 



6 ' 12 

Co 6„ + Cl 6„_i + C2 6„_2 + C3 6„_3 



C4 6„_4 



(n + l)(n + 2) 



20 



with n > 4. 



(51) 
(52) 
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From ([35|l and ([37|) we find (we use C2 = 1) 

60 = v{ax), bi = ifi'ia^). (53) 

Thus, using the coefficients &o and 61, one can define the wave function and its derivative at the point a^- Defining bo 
and &i and using ([5^ . one can calculate all &„ needed. 

Analysis has shown that this way is more efficient (in comparison with the previous one) in the calculation of the 
decreasing component of the wave function. Therefore, using the escape point atp as the starting point, we shall use 
this way to calculate the second partial solution, which is decreasing in the barrier region, and we select 

bo = 1, 61 = -0.1. (54) 

3.2.3 Determination of the first partial solution in the whole region, increasing in the region of tunneling 

Picking values for the wave function and its derivative at a = and calculating them in the vicinity of that point, we 
shall calculate them in the whole range of a, using the Numerov method with a constant step. We see that such an 
approach allows one to obtain a convergent and stable partial solution, which is increasing in the tunneling region. In 
order to analyze to what extent this approach gives convergent (stable) solutions, for a comparison we use |23j with 
data presented for the modulus of the wave function (see (9) p. 5; we use the parameters of the potential in (|13p). 

The first partial solution for the wave function and its derivative, obtained by our calculations in the region of 
the escape point, is shown in Fig. [21 Comparing them with results in Ref. [23 1 we see that in our approach the wave 
function and its derivative look more continuous and have no divergences (this becomes more apparent after a detailed 
analysis of the selected regions). In the figures one can see that starting from the point atp the wave function decreases 
and its derivative increases with increasing a. One can see that the derivative looks significantly larger than the wave 
function. 



5,00E+261 




Parameter a 

Fig. 2. The first partial solution for the wave function and its derivative (decreasing in the tunneling region): curve 1, blue, is 
for the wave function; curve 2, green, for the derivative of this wave function. 

Now let us consider the behavior of the wave function more carefully. In Fig. [3] it is shown how it changes in the 
tunneling region and in the vicinity of the escape point. From the figures one can see that the wave function satisfies 
the rules of behavior of the wave functions inside sub-barrier and above the barrier region. Starting from very small 
a, the wave function increases monotonously (without any oscillations) with increasing a, which corresponds to the 
tunneling region (this becomes more obvious especially in a logarithmic presentation of the wave function; see the left 
figure). Further, oscillations have appeared in the wave function, which can be possible only inside the above barrier 
region (a smooth transition is shown in the right figure). The boundary of such a transition of the wave function must 
be the escape point atp (our calculations give atp = 17.320508, which coincides with Ref. [23]). 

In Fig. |4] it is shown how the wave function and its derivative vary at large values of a, i. e. including a region 
which may be called an asymptotic one. We obtain smooth continuous solutions running up to a = 100 (in Ref. |23| 
the maximum value is a = 30). An interesting peculiarity of the obtained solutions is an absolutely uniform straight 
increase of the maxima of the derivative of the wave function and a smooth decrease of the wave function with 
increasing a, which must point to the specific character of the expansion of the Universe in this quantum-mechanical 
approach. 

Analyzing all the obtained figures, we find the following interesting property. 
With increasing a, the period of oscillations decreases uniformly both for the wave function and its derivative in the 
above barrier region! 
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Fig. 3. The behavior of the wave function in the selected regions: (a) is for the wave function in the tunneling region (In {(fii{a)) 
is in vertical axis); (b) for the wave function in the region of the escape point. 




— This explains why for sufficiently small increases of a it becomes significantly more difficult to calculate the con- 
vergent continuous solutions for the wave function and its derivative (where it is needed to significantly decrease 
the step)! 

— Instead of the Numerov method with constant step, it can be more efficient to use another method of continuation 
of the solution, in which the step decreases uniformly with increasing of a. 

— This information can be interesting for the choice or construction of new functions, by which the total wave function 
in the interesting region can be expanded with the highest efficiency (for example, this is why searching for solutions 
for the total wave function as an expansion in plane waves or spherical Bessel functions cannot give a sufficiently 
stable convergent result (with a small increase of the selected range of a)). 



3.2.4 Determination of the second partial solution on the whole region, decreasing in the tunneling region 

We find the second partial solution similarly to the first one, but with a different starting point. First, setting the 
value for the wave function and its derivative at the escape point atp, we find them in the vicinity of that point by 
the method from sec. 13.2.21 Further, using the Numerov method with a constant step, we find solutions in the whole 
region of a. 

In Fig. [5] the wave function and its derivative in the region with the right boundary larger than point atp are shown 
(parameters of the potential are from ([13])). One can see from this figure a strong decrease of the amplitudes of the 
wave function and its derivative with increasing a, which implies that the solutions obtained are decreasing in the 
tunneling region. 

Let us analyze the behavior of the wave function more carefully. In Fig. [6] its shape in the region of tunneling of 
the wave and close to the escape point is shown. In the first figure one can see its monotonous logarithmic decrease 
(without oscillations) in the tunneling region with increasing a. In the second figure the behavior of the wave function 
in the region of the escape point is shown, where the first oscillations appear. In the third figure one can see how a 
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12 14 15 18 20 



Parameter a 



Fig. 5. The second partial solution for the wave function and its derivative, decreasing in the tunneling region: curve 1, blue, 
is for the wave function; curve 2, green, for the derivative of this wave function. 



smooth decrease of the wave function is transformed little by little into the first oscillation in the small neighborhood 
of point atp, and the boundary of such a transformation must be the escape point atp (as in the first solution, here 
the rules of behavior for the wave function in the sub-barrier and above barrier regions are shown too). As we see, the 
proposed approach gives smooth pictures without divergences both for the wave function and its derivative. 




Parameter a Parameter a Parameter a 

Fig. 6. The behavior of the wave function in the selected region: (a) the wave function in the tunneling region (for the vertical 
axis we use In {1^2(0.))); (b) the wave function in the region of the escape point; (c) the smooth transformation of the decrease 
of the wave function into the first oscillation close to the point atp. 



Analyzing all pictures, we point out the following properties. 

— The principally different behavior of the two found partial solutions for the wave function (and its derivative) in 
the tunneling region ensures their linear independence. This allows one to construct a stable general solution for 
the wave function and for its derivative on their basis. 

— A property of the uniform decrease of the oscillation period in the above barrier region for the wave function and its 
derivative with increasing a, found for the first solution (increasing in the tunneling region), exists for the second 
solution (decreasing in the tunneling region) as well. 

— As this property is satisfied for each linearly independent partial solution, it must be satisfied for the total wave 
function and its derivative as well. From this one can obtain new information on the .specific character of the 
dynamics of the Universe expansion (and to use the oscillation period as a new characteristic for the description 
(estimation) of this dynamics). 

— Two partial solutions for the wave function and its derivative and their analysis are obtained at zero energy. 



3.3 Determination of the total wave function in the whole semiaxis of a 



After obtaining two linearly independent partial solutions (pi{a) and ^'2(0) for the wave function at the potential l|12[) . 
we write the general solution as 

ipia)^N-{Ciipi{a)+C2Ma-)), (55) 
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where Ci and C'2 are arbitrary complex constants, is a normalization factor (we separate it explicitly from the 
constants Ci and C2, because these constants and the factor N are found from different conditions, and they have 
different effects on the penetrabihty coefficient and other characteristics). 

One can find the unknown constants Ci and C2 from the tunnehng boundary condition pointed out above: the 
function f{a) must define only an outgoing wave at the escape point atp. Here, the function ip{a) and its derivative 
must be equal to ([32]) at the point atp (we shall use the factor N for (p{a) in consistence with the factor N in ([32]) ): 



<fiatp) = NW{^ = 0), 



dip{a) 



da 



N 



da 



N\Vi 



,1/3 dm 



d^ 



(56) 



Taking into account eq. (|55p . we find (a prime is for a derivative with respect to a) 



Ci 



, C2 - 



(57) 



The total wave function for the potential with the parameters from (|f 3[) . calculated by our approach, is shown in 
Fig. El From Fig.[5]one can analyze the behavior of the modulus of this wave function and from Fig.[Tn] — the behavior 




Parameter a Parameter a 

Fig. 7. The total wave function: curve 1, blue, is for the real part; curve 2, green, for the imaginary part. 



of its phase. 




Parameter a Parameter a Parameter a 



Fig. 8. The modulus of the total wave function: (a) in the tunneling region the modulus decreases uniformly as a function 
of a; (b) in the escape region, the minima and maxima appear in the wave function, which is connected with the existence of 
oscillations in it, but the modulus is never equal to zero in the whole range of a (this demonstrates the existence of a non-zero 
flux); (c) near the escape point atp we see that the modulus is changed minimally under increasing of a (this demonstrates the 
fulfillment of the definition introduced above for the wave at point atp). 
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Parameter a Parameter a Parameter a 



Fig. 9. The phase of the total wave function: (a) in the barrier region the phase is constant, outside the barrier oscillations 
appear in the limits from — tt to vr; (b) the smooth appearance of the first oscillations of the phase in the neighborhood of the 
escape point atp is shown; (c) the oscillation period of the phase decreases uniformly with increasing a. 

Table 1. Probability Ptun in dependence on the upper limit amax and decomposition of the whole range of a. 





Decomposition 


10000 


50000 


100000 


200000 


1000000 


2000000 


20 


0.0428858066 


0.0420408533 


0.0418943250 


0.0419056013 


0.0418471293 


0.0418398319 


30 


0.0230343856 


0.0229506126 


0.0230091568 


0.0229690161 


0.0229369614 


0.0229329589 


40 


0.0197653508 


0.0197616584 


0.0195350357 


0.0195796586 


0.0195208513 


0.0195095816 


50 




0.0183045342 


0.0180612516 


0.0180298553 


0.0180046869 


0.0180015439 


70 






0.0167357966 


0.0165972606 


0.0165979553 


0.0165950583 


100 








0.0156822095 


0.0156892697 


0.0156865629 


150 










0.015056842 


0.0150574352 


200 












0.0147582032 



4 Estimation of the barrier penetrability 

For the determination of the barrier penetrability in standard quantum mechanics for the wave escaping from this 
barrier outside we must know the incident flow inside the internal region. This flow must be constructed on the basis 
of the wave already tunneling under the barrier starting from point a = to the right. This method requires an 
accurate analysis and some questions may appear. Therefore, with the purpose to estimate the penetrability of the 
barrier studied, we restrict ourselves only to two other characteristics (which are analyzed often in the development 
of quantum cosmological models): we define the probability Ptun of the appearance of the particle studied inside the 
tunneling region and the probability Pcxt of the absence of it as the corresponding ratios of the probability of the wave 
functions, defined on the tunneling region and on the external region, to the probability of the wave function, defined 
over the whole range of its definition: 

Qmax atp 

I Ma)\^da J \^{a)\-'da 

Pext = , -Ptun = (58) 

Y\^{a)\'da Y\ip{a)\^da 



where amax is the upper limit of the range of a where we consider the total wave function ^p{a). Note that we use 
the stationary wave function Lp{a) at zero energy in contrast to [23j . In order to keep the computer calculations as 
accurate as possible, we calculate Pgxt thus: 

Pe,t = l-Z'tun. (59) 

The results of our calculations of Ptun are presented in Table [1] One can see that increasing of the decomposition 
in the selected region gives a stable convergent value for Ptun- This demonstrates the convergence of the proposed 
method in the calculation of the wave function. From the data we select the following limit: 

Ptun = 0.014, Pe,t = 1 - 0.014 = 9.986, (60) 

which differs signiflcantly from the results of |23j . 



16 



Sergei P. Maydanyuk: Wave function of the Universe in the early stage of its evolution 



4.1 Comparison with semiclassical calculations 



Now we shall estimate the penetrability in the semiclassical approximation. First, let us consider the wave function 



in the case E ^ {]. Here, we have internal and external turning points, which we denote a^^ 



semiclassical wave function in the region to the right from the external turning point a^^ 
the barrier outside having the form (we define it in the second approximation) 



<Pout(a) = N ■ f2{a) ■ exp 



(cxt) 




be a wave outgoing from 



(61) 



where /2(a) — \/ \/\p (a)|, p{a) — \l E — V^(a) is the complex momentum, and A'^ is a normalization factor. Then 
according to the rules of correspondence (50,2) and (47,5) in [16^ between wave functions in different regions close to 
the turning points, we have obtained the wave function (ptun(i) in the barrier region and the wave function (pint(i) in 
the left region before the barrier: 



'/'tun(a) ^ N ■ f2{a) ■ exp 



(cxt) 



p da 



TV •/2(a)- exp ^ 



(ext) 



p da 



(int) 
tp 
(int) 
tp 



p da 



(int) 



^(cxt) ^Ont) 

(pi„t(a) = 2 iV • /2(a) • exp ( J \p\ da j ■ coal J pda-^\ 

\(int) ^ a 



(62) 



In Fig. [To] one can see the modulus of the wave functions calculated in the first and second semiclassical approxi- 
mations (at E = 0) located together with the modulus of the wave function obtained by the direct quantum method 
proposed above. Here, we calculate the wave function in the second semiclassical approximation inside the barrier 
region by the first formula in (|62[) and in the external region by the second formula in (|62[) . The wave function in 
the first semiclassical approximation we calculate by the same formulas and use /2 = 1. We normalize both solutions 
found separately on the wave function obtained by the direct approach, using the factor N. From these figures we 




Fig. 10. The modulus of the wave function calculated by the direct quantum approach and wave functions calculated in the 
first and second semiclassical approximations: (a) in the consideration inside the whole barrier region, we obtain a complete 
and very precise coincidence between the modulus of the wave function calculated by the direct approach and both wave 
functions calculated by two semiclassical approximations (except for a region close to zero for the solution in the second WKB 
approximation, which singularity is not visible at such scale); (b) close to zero the wave function in the first WKB approximation 
describes more accurately the modulus of the wave function in the direct approach in comparison with the wave function in the 
second WKB approximation (at a = 0, the wave function in the second WKB approximation has a singularity, while we obtain 
finite values for the wave function in the direct approach and the wave function in the first WKB approximation); (c) near the 
escape point atp the modulus of the wave function in the second semiclassical approximation coincides with the maxima of the 
wave function in the direct approach. 



conclude to the following. 

— Inside the whole barrier region the modulus of the wave function in the direct approach is approximated very 
precisely and accurately by the wave functions in the first and second approximation. Close to zero it is described 
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better by the wave function in the first approximation: for the wave functions in the second approximation we have 
a singularity while for the wave function in the direct approach and for the wave function in the first approximation 
we have obtained finite values at point a = 0. In a detailed study of three solutions we observe their values to 
be different (with a similar tendency), which can be shown in further calculations of the penetrability and similar 
characteristics. By this, one can explain different possible (estimated) values of the penetrability at the very precise 
coincidence between the exact and semiclassical solutions presented. 

— Inside the external region the wave function in the first approximation can be considered as a "plane" wave (with 
constant modulus and changed phase). We see that it describes less accurately the exact solution while the wave 
function in the second approximation coincides very accurately with maxima of the wave function in the direct 
approach, and it can be considered as a proper solution. We observe the presence of oscillations of the modulus 
of the wave function in the direct approach which decrease with increasing of a (and which are not present in 
semiclassical solutions). Such oscillations may indicate the following supposition: the outgoing wave in the external 
region close to the turning point cannot be considered as a usual "plane " wave, and oscillations include information 
of the interaction between the studied outgoing wave and the barrier at not too large distances. One can connect 
them with the phenomenon of the barrier non-locality: the barrier has influence on the tunneling and further 
propagating wave as a unified object, and in estimating of the penetrability we consider the barrier as the potential 
in the whole region. The rules of the correspondence between two semiclassical solutions of the wave function 
in different regions can break down (or change) the non-locality of the barrier near the turning point, and we 
obtain a greater local influence of the barrier on the semiclassical solutions of the wave function in different regions 
separately. 

— Close to the turning point the wave function in the direct approach connects both semiclassical solutions for the 
wave function in tunneling and external regions. Such a connection can be considered as a small shift of the wave 
function along the a axis (perhaps this has been found for the flrst time in cosmological models). 



Now on the basis of ([61]) and ([62|) we find the known formula for the penetrability coefficient in the semiclassical 
approximation: 



(oxt) 



^'penctrabmty=exp -2 / \p\ da ] . (63) 



(i„t) 



Taking into account only the first semiclassical approximation and comparing the first expression in (|62p with (|63p 
(at the needed normalization of N), we obtain a new formula: 



pWKB, (1) 
penetrability 



f (int)s 



(64) 



/ |^(a)Pda 





Calculations in the semiclassical approximation gives (decomposition of the barrier region is in 100000 intervals): 

= 0.052, Pp^eLbiiity = = 7.0 • IQ-^^^ (65) 

Now we see that the direct approach gives the i^t'u"'^* coefficient 4-5 times smaller than P^^^. This points to the 
existence of more essential tunneling processes inside the barrier. One may also suppose that taking the external region 
into account we can change the value of the penetrability coefficient ^'jj^j^^iabiiity found. 



5 Conclusions and perspectives 

In this paper a new method of (non-scmiclassical) calculation of the wave function of the Universe in the quantum 
cosmological model of a Universe of closed type, in the framework of the Friedmann-Robertson-Walker metric, is 
presented. We note the following. 

— A method for the calculation of two partial solutions for the wave function and its derivative is constructed: 

— at first, the values for the wave function and its derivative are defined at a starting point (for the partial solution 
which increases in the barrier region, we use the starting point a = 0, and for the second partial solution, which 
decreases in the barrier region, we select the starting point to be the escape point atp); 

— using the algorithms in Sees. 13.2.11 and 13.2.21 both partial solutions are found independently in the vicinity of 
the starting point; 
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— using the Numerov method with a constant step, both partial solutions are found in the whole studied range 
of a. 

In this way, we obtain a well-convergent picture for the wave function and its derivative (for the potential from 
Ref. [23], the convergence of the solutions is demonstrated up to a < 100). 

— In order to describe the process of tunneling as accurately as a possible, to construct the total wave function 
on the basis of its two partial solutions unambiguously, we use the tunneling boundary condition from papers of 
A. Vilenkin, at the escape point atp- the total wave function must represent only the outgoing wave. 

— We introduce two definitions of this wave in the (right) neighborhood of the boundary atp as follows. 

— The wave is such a linear combination of two partial solutions of the wave function that a change the modulus 
p of this wave function is the closest to constant under variation of a (strict definition of the wave; see (HH)) 

— The wave is such a linear combination of two partial solutions of the wave function that the modulus p changes 
minimally under variation of a (weak definition of the wave; see (|22p ). 

— Such definitions of the wave allow us to construct sufficiently stable and convergent solutions for the total wave 
function and its derivative (for the potential from [23| the solutions are calculated at a < 100). 

— Analyzing the wave function, the following properties have been shown. 

— We observe oscillations of the modulus of the wave function in the external region starting from the turning 
point atp which decrease with increasing of a, and these are not shown in semiclassical calculations (perhaps, 
they have been found for the first time). An interesting idea is in the physical interpretation: their maxima an 
minima point to values of a where the modulus of the wave function of the Universe is maximal or minimal, 
and one may suppose that such an effect shows specific peculiarities of space-time of the Universe in smaller 
distances while it disappears at large ones, tending to the semiclassical limit. 

— With increasing a, the period of oscillations decreases uniformly smoothly both for the wave function, and for 
its derivative in the above barrier region. This peculiarity is fulfilled for each partial solution, and, therefore, it 
must appear in the total wave function and its derivative as well. From this one can find new information on 
the specific character of dynamics of the expansion of the Universe (and one can use the oscillation period as 
a certain characteristic for the estimation of the dynamics of such an expansion). 

• This explains why for a sufficiently small increase of a it becomes significantly more difficult to calculate 
the convergent continuous solutions for the wave function and its derivative (where it is required to reduce 
the step significantly)! 

• This information may be interesting for the choice or construction of new functions, by which the total 
wave function in the region of interest can be expanded with the highest efficiency (for example, this is 
why searching for solutions for the total wave function as an expansion in plane waves or spherical Bessel 
functions cannot give a sufficiently stable convergent result (with a small increase of the selected range of 
a))- 

• In order to calculate the wave function and its derivative with higher efficiency, we propose, instead of the 
Numerov method with a the constant step, to use another method of continuation of the solution, where 
the step decreases uniformly with increasing a. 



Acknowledgements 

The author appreciates useful discussions with Dr. Volodymyr Uleshchenko concerning the application of quantum non- 
semiclassical approaches in quantum cosmology and determination of wave functions of the Universe inside the tunneling 
region, and he is graceful to Dr. Yevgeny Kats for his assistance in preparing the paper, and Dr. Mariam Bouhmadi Lopez for 
useful comments concerning different aspect of quantum cosmology of asymptotically de Sitter Universes. 



References 

1. A. Vilenkin, Predictions from quantum cosmology, Proceedings: String Gravity and Physics at the Planck Energy Scale 
(International School of Astrophysics (NATO Advanced Study Institute), Erice, Italy, 8-19 Sep 1995; NATO ASI Series 
C, Mathematical and Physical Sciences, 476; Edited by N. Sanchez and A. Zichichi; Boston, MA, Kluwer, 1996, 544 p.), 
p. 345-367, gr-qc/9507018 

2. B. S. DeWitt, Quantum theory of gravity. I. The canonical theory, Phys. Rev. 160 (5), (1967) 1113-1148. 

3. J. A. Wheeler, Batelle Rencontres (Benjamin, New York, 1968). 

4. A. Vilenkin, Creation of universes from nothing, Phys. Lett. B 117 (1-2), (1982) 25-28. 

5. J. B. Hartle and S. W. Hawking, Wave function of the Universe, Phys. Rev. D 28 (12), 2960-2975 (1983). 

6. A. D. Linde, Quantum creation of the inflationary Universe, Lett. Nuov. Cim. 39 (2), (1984) 401-405. 



Sergei P. Maydanyuk: Wave function of the Universe in the early stage of its evolution 



19 



7. Ya. B. Zel'dovich and A. A. Starobinsky, Quantum creation of a universe in a nontrivial topology, Sov. Astron. Lett. 10 
(3), (1984) 135. 

8. V. A. Rubakov, Quantum mechanics in the tunneling universe, Phys. Lett. B 148 (4-5), (1984) 280-286. 

9. A. Vilenkin, Quantum creation of universes, Phys. Rev. D 30 (2), (1984) 509-511. 

10. A. Vilenkin, Boundary conditions in quantum cosmology, Phys. Rev. D 33 (12), (1986) 3560-3569. 

11. P. L Fomin, Dokl. Akad. Nauk Ukr. SSR 9A, (1975) 831. 

12. D. Atkatz and H. Pagels, Origin of the Universe as a quantum tunneling effect, Phys. Rev. D 25 (8), (1982) 2065-2073. 

13. A. Vilenkin, Approaches to quantum cosmology, Phys. Rev. D 50 (12), (1994) 2581-2594, gr-qc/9403010 

14. V. A. Rubakov, Quantum cosmology, Proceedings: Structure formation in the Universe (Edited by R. G. Crittenden and 
N. G. Turok, Kluwer, 1999), p. 63-74, gr-qc/9910025 

15. J. Ambjorn, J. Jurkiewicz and R. Loll, Semiclassical universe from first principles, Phys. Lett. B 607, (2005) 205-213. 

16. L. D. Landau and E. M. Lifshitz Quantum mechanics, course of Theoretical Physics, Vol. 3 (Nauka, Mockva, 1989), p. 768 

- [in Russian; eng. variant: Oxford, Uk, Pergamon, 1982]. 

17. S. Weinberg, Gravitatsiya i kosmologiya: printsipi i prilozheniya obschei teorii otnositel'nosti (Mir, Moskva, 1975), 696 p.. 

- [in Russian; eng. variant: S. Weinberg, Gravitation and cosmology: principles and applications of the General theory of 
relativity, MIT, John Wiley and Sons, New York - London - Sydney - Toronto, 1972]. 

18. V. A. Rubakov, Introduction to cosmology (RTN Winter School of Strings, Supergravity and Gauge Theories, January 31 

- February 4 2005, SISSA, Trieste, Italy), 58 p., - available in (pos.sissa.it | 

19. A. Linde, Particle physics and inflationary cosmology (Harwood, Chur, Switzerland, 1990), 362 pp.; Contemporary 
Concepts in Physics, v. 5, hep-th/0503203 

20. M. Trodden and S. M. Carroll, TASI Lectures: Introduction ot cosmology. Lectures at the Theoretical Advanced Study 
Institutes in elementary particle physics (TASI-2003, Recent Trends in String Theory, University of Colorado at Boulder, 
1-27 Jun 2003; Edited by J. M. Maldacena; Hackensack, World Scientific, 2005; 548 pp.), 82 pp., |astro-ph/0401547[ 



21. R. H. Brandenberger, Inflationary cosmology: progress and problems. Lectures at the International School on Cosmology 
(Kish Island, Iran, Jan. 22 - Feb. 4 1999; Proceedings: Kluwer, Dordrecht, 2000), p. 48, hep-ph/9910410 

22. D. Levkov, C. Rebbi and V. A. Rubakov, Tunneling tn quantum cosmology: numerical study of particle creation, Phys. 
Rev. D 66 (8), (2002) 083516, gr-qc/ 0206028 , 

23. J. Acacio de Barros, E. V. Correa Silva, G. A. Monerat, G. Oliveira-Neto, L. G. Ferreira Filho and P. Romildo Jr., Tunneling 
probability for the birth of an asymptotically de Sitter universe, Phys. Rev. D 75, (2007) 104004, 'gr-qc/0612031 , 

24. V. V. Kuzmichev, Evolution of the quantum Friedmann Universe featuring radiation, Phys. Atom. Nucl. 62, (1999) 708-714; 
Yad. Fiz. 62, (1999) 758-764, gr-qc/0002029 

25. V. E. Kuzmichev and V. V. Kuzmichev, Properties of the quantum universe in quasistationary states and cosmological 



puzzles, Eur. Phys. Journ. C 23, (2002) 337-348, astro-ph/01 11438 



26. Spravochnik po spetsialnim funktsiyam s formulami, grafikami i matematicheskimi tablitsami. Pod redaktsiei M. Abramow- 
itza i I. A. Stegan (Nauka, Moskva, 1979), 832 p.. - [in Russian; eng. variant: Handbook of mathematical functions with 
formulas, graphs and mathematical tables. Edited by M. Abramowitz and I. A. Stegan, National bureau of standards. 
Applied math, series - 55, 1964]. 



